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ABSTRACT 

Based on a 21.5 ks Chandra observation of A2556, we identify an edge on 
the surface brightness profile (SBP) at about 160/if/ kpc northeast of the cluster 
center, and it corresponds to a shock front whose Mach number Ai is calculated 
to be 1.25;to;o3- prominent substructure, such as sub-cluster, is found in either 
(-H I optical or X-ray band that can be associated with the edge, suggesting that the 

Q-[ conventional super-sonic motion mechanism may not work in this case. As an 

alternative solution, we propose that the nonlinear steepening of acoustic wave, 
which is induced by the turbulence of the ICM at the core of the cluster, can 
be used to explain the origin of the shock front. Although nonlinear steepening 
weak shock is expected to occur frequently in clusters, why it is rarely observed 
still remains a question that requires further investigation, including both deeper 
X-ray observation and extensive theoretical studies. 

Subject headings: shock waves — galaxies: cluster: individual (Abell 2556) — 
(Galaxies:) intergalactic medium — X-ray: galaxies: clusters 



1. INTRODUCTION 

Radiative loss via X-ray emission in galaxy clusters has long been supposed to be the 
origin of the substantial gas inflow, namely the "cooling flow" (see Fabian 1994 for a review). 
However, X-ray spectra from ASCA and XMM-Newton observations fail to show line emission 
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from ions with intermediate or low temperatures, which imphes that the coohng rate is only 
at most 20% of the previously assumed value (e.g., Kaastra et al. 2001; Peterson et al. 2001; 
Xu et al. 2002; Peterson et al. 2003; Tamura et al. 2003). Chandra observations have 
also confirmed such small cooling rates (McNamara et al. 2000). There must be a heating 
source in the cluster core region. The most popular candidate of the heating source is the 
AGN at cluster centers. The AGN dissipates energy via bubbles, jets, sound waves and weak 
shocks evolved from sound waves (Fabian & Sanders 2009). In Perseus and Virgo Clusters, 
previous works have observed the sound waves and weak shocks generated by activities of 
central AGNs (Fabian et al. 2003; Ruszkowski et al. 2004; Forman et al. 2005). Besides the 
AGN activity, some other heating mechanisms have also been concerned in the last decade. 
Especially since the launch of Chandra, shock fronts have been revealed clearly in the ICM 
of galaxy clusters through X-ray imaging-spectroscopic studies, which can efficiently convert 
kinetic energy of the gas into thermal energy. Such heating process may be associated with 
the feedback mechanism that balances the radiative cooling in galaxy clusters. 

This work would study the weak shock front in the galaxy cluster A2556 by using the 
Chandra observation. A2556 (RA=23hl3m03.3s, Dec=-21d37m40s, J2000.0; richness 1) is 
located at z ~ 0.0871 in the Aquarius supercluster (Batuski et al. 1999), an overdense region 
of the universe. It is moderately luminous in the X-ray band (Lx ^ 2 x 10^^ erg s~^; Ebeling 
et al. 1996), possessing a velocity dispersion ~ 872 km s~^ (White et al. 1997) and a virial 
mass ^ 2.5 ± 0.1 x 10^^ (Reimers et al. 1996). 

We describe the data reduction, imaging and spectroscopy analysis in §2; we discuss 
and summarize our results in §3 and §4, respectively. Throughout this paper, we adopt the 
cosmological parameters Hq = 71 km s~^ Mpc^^, = 0.27 and I^a = 0.73, so that 1' 
corresponds to 96.7/if/ kpc at such distance. We utilize the solar abundance standards of 
Grevesse and Sauval (1998), where the iron abundance relative to hydrogen is 3.16 x 10~^ in 
number. Unless stated otherwise, all quoted errors are derived at the 68% confidence level. 

2. DATA ANALYSIS 

2.1. Observation and Data Reduction 

The Chandra observation of A2556 was carried out on October 5, 2001 (ObsID 2226) 
for a total exposure of 21.5 ks with CCDs 3, 5, 6, 7 and 8 of the Advanced CCD Imaging 
Spectrometer (ACIS) in operation. Events were collected with a frame time of 3.2 s and 
telemetered in the VeryFaint mode. The focal plane temperature was set to —120 °C. The 
center of A2556 was positioned nearly on the ACIS-S3 chip with about 0.578' offset. Most of 
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the X-ray emissions of A2556 are covered by the S3 chip. Therefore, this work focuses on the 
data drawn from the S3 chip. The Chandra data analysis package CIAO software (version 
3.4) is used to process the data from the S3 chip. We keep events with ASCA grades 0, 2, 3, 
4 and 6, and remove all the bad pixels, bad columns, columns adjacent to bad columns and 
node boundaries. No background flare is found during the live operating time. 



2.2. X-ray Surface Brightness 

We show the raw Chandra ACIS-S3 image of A2556 in 0.7-7.0 keV band in Figure [1^, 
and have all the point sources excluded, which could be detected at the confidence level 
of 3(7 by the CIAO tool celldetect. The black cross in the center of Figure [1^ indicates 
the X-ray peak of A2556. The core region of the A2556 shows no complex irregularity and 
should be composed of a single nucleus; the outer region shows a slight elongation toward the 
northwest. The X-ray morphology roughly shows a generally relax appearance, except that 
there exhibits an arm-like structure with an edge at approximately IGOhj^ kpc northeast of 
the cluster center, and the edge spans about 85° azimuthally from southeast to northeast. 
We highlight the edge span by two white lines in Figure [1^. Figure [TJj shows the SDSS B- 
band image of A2556 with the Chandra 0.7-7.0 keV X-ray intensity contours overlaid. The 
X-ray peak consists with the centroid of the cD galaxy 2MASX J23130142-2138039 within 
1". The cD galaxy shows an elongation in the same direction as the X-ray morphology, which 
may due to the nearby A2554. 

To get more insight into the edge of this arm-like structure, we extract the exposure- 
corrected surface brightness profile (SBP) from a series of narrow elliptical annuli, which 
are parallel to the east (E) region show in the Figure [Dd. In order to get a comparison, 
we also extract SBPs of other three directions in Figure [T)d, namely south (S), west (W) 
and north (N) regions. The energy band for these profiles is restricted to 0.7-7.0 keV. We 
show such extracted SBPs in Figure [21 In region E, a substantial discontinuity of the SBP 
could be observed at approximately IGOh^^ kpc, and other three regions do not present such 
characteristic. 

To quantitatively describe the discontinuity of the SBP in region E, we fit observed SBP 
with a modeled gas emissivity distribution, which corresponds to a function as 



where r represents the radius of the projected SBP; the gas emissivity distribution, e{r), is 
broken at the truncation radius, -Rcut, and described with two power-law components; ei, 
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£2, tti and 0^2 are 4 fitting parameters. The solid line in Figure |2] shows the best-fit SBP 
obtained with Eq. [TJ The best-fit -Rcut = 158.26lQ;42/if/ kpc corresponding to the edge of 
the arm-like structure, ei and 62 shall represent the plasma emissivities inside and outside of 
the edge, respectively. The best-fit emissivity jump ei/e2 = 1.901q;q4. Since the gas density 
is related to the emissivity by e ~ p^, the density jump across the edge is calculated to be 
P1/P2 = {er/e^f' = l-STt^i 



2.3. Spectral Analysis 

In the Chandra spectral analysis, there are two usual ways to obtain a background 
spectrum, including cosmic, instrumental and non X-ray particle components. Firstly, the 
observated dataset can be directly used to extract the background spectrum from a local 
and uncontaminated region on Chandra chips, which must be located far away enough from 
the source. Secondly, one can draw the background spectrum from the Chandra blank field 
datasets in the same detector region as where the source spectrum is extracted. Since the 
core radius of A2556 is 350/?,^/ kpc (White et al. 1997) and it should cover most of ACIS-S3 
CCD, we have to utilize the Chandra blank-sky template as the background for A2556. The 
template is tailored to match the actual pointing and the background spectrum is extracted 
and processed identically to the source spectrum. Then, the background spectrum is rescaled 
by normalizing its high energy end to the corresponding observed spectrum. Corresponding 
spectral redistribution matrix files (RMFs) and auxiliary response files (ARFs) are created 
with CIAO tools, mkwarf and mkacisrmf , respectively. All spectra are rebinned to give 
a minimum of 20 raw counts per spectral bin, which allows statistics to be applied. 
Since the contribution of the hard spectral component is expected to be rather weak, and 
also to minimize effects of the instrumental background at higher energies and calibration 
uncertainties at lower energies, the spectrum fitting is restricted to 0.7-7.0 keV . 



2.3.1. Temperature Across the Edge 

In order to derive the temperature profile across the SBP edge, we extract the spectra 
from 4 elliptical annuli in region E of Figure [T]d, so that the SBP edge lies exactly at 
the boundaries of certain elliptical annulus, and the elliptical annulus covers the azimuthal 
angle, where the edge is most prominent. We use the XSPEC 12.4.0 package to perform the 
deprojected spectral analyses. In spectral fitting process, the gas emission is modeled with 
an absorbed APEC component and the PROJCT model, which performs a 3-dimensional 
to 2-dimensional projection of prolate ellipsoidal shells onto elliptical annuli and evaluates 
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the influence of outer ellipsoidal shells onto inner ones. The absorption column density iVn 
is fixed at the Galactic value A^h = 2.07 x 10^° cm~^ (Dickey & Lockman 1990). Best- 
fit temperature values are shown in Figure [3^. Temperature behind the edge of arm-like 
structure is much higher than those of ambient regions. The calculated temperature jump 
across the edge is T1/T2 = 1.931!q;72- Both gas density and temperature jumps across the 
edge of the arm-like structure indicate that the edge is a weak shock front. 

The Mach number, Ai, of the weak shock front can be determined by employing the 
Rankine-Hugoniot shock relation (Landau & Lifshitz 1959) 

p, ^ (1 + 7)M^ 

P2 2 + (7 - 1)M2 ^ ' 

and 

Ti ^ (27M^ - (7 - 1)) (2 + (7 - 1)M') 

T2 (7 + 1)2M2 ' ^ ^ 

where Ti and T2 denote the temperature behind (the "post-shock" ; region 1 in Figure [TJd) and 
ahead (the "pre-shock" ; region 2 in Figure [T]d) of the shock, respectively; pi and p2 denote 
the density behind and ahead of the shock, respectively. For the intracluster gas, 7 = 5/3, 
which is the adiabatic index for monatomic gas. For weak shock fronts, the accuracy of Ai 
derived from the density discontinuity is better than that from temperature discontinuity 
(Markevitch & Vikhlinin 2007). Therefore, with the density jump pi/ p2 = 1.37too3 -^l- 
[21 we derive the Mach number Ai = 1.25l;g o3. According to Eq. [3], a temperature jump of 
T1/T2 = 1.241qq^ is indicated, which confines the spectroscopic value 1.931^072 within errors. 

We crosscheck the detected temperature jump by examining the 2-dimensional tempera- 
ture distribution of the gas in A2556. In Figure [3]d, the projected temperature map of A2556 
(Gu et al. 2009) exhibits a substantial high temperature region exactly behind the edge of 
the arm-like structure, which is coincident with the temperature profile shown in Figure [3^. 



2.3.2. Classification of Weak Cool Core Cluster 

To test whether a cooling flow model fits the spectrum of A2556, we apply four types of 
models to fit the spectrum in a circle region with radius of 237/1^/ kpc, which center on the 
X-ray peak. The circle region covers nearly the whole Chandra ACIS-S3 CCD. Firstly, we 
use the WABSxAPEC model (hereafter ITA) and WABSxMEKAL model (hereafter ITM) 
to describe the X-ray emission only by a single isothermal plasma. In these two models, 
temperatures and elemental abundances are allowed to vary. Secondly, in order to feature 
a cooling flow, we employ a single temperature MEKAL model plus a multi-phase cooling 



- 6 - 



flow component MCKFLOW, where their abundances are tied together and allowed to vary; 
the temperature of MEKAL and the higher gas temperature of MCKFLOW are also tied 
together and allowed to vary. We consider two cases for this model: the lower bound on the 
gas temperature in the MCKFLOW is fixed at 0.1 keV, which means the cooling material is 
fully cooled below the X-ray emitting temperature (hereafter denote FC for fully cooling), 
and this component is set as a free parameter (hereafter denote VC for variable cooling). 

Best-fit spectral parameters of four models with 90% confidence level are presented in 
Table 1. As a general comment, all models do a reasonably good job at fitting the spectrum 
with reduced values around 1.1. The ITA and the ITM give nearly the same spectral 
parameters. The cooling rates for the FC and the VC are 19.99l}|];}^ Mq yr'^ and 20.13tl%f 
Mq yr~^, respectively. Therefore, A2556 could be classified as a weak cool core cluster. This 
conclusion consists with the cooling time within 1.5' of A2556, 6.5 Gyr (Gu et al. 2009), 
according to which A2556 could also be classified as a weak cool core cluster (Mittal et al. 
2009). 



3. DISCUSSION 

In §2, we present robust evidences for the existence of the weak shock front located at 
approximately IGO/i^/ kpc northeast of the cluster center. As in the Perseus, the Hydra A, 
and several other clusters, the AGN activity might be the cause of the weak shock in A2556. 
However, partly due to the poor photons statistics, we fail to find robust evidence of AGN 
bubbles and disturbances on cluster scales. In this section, we shall discuss the other two 
possible origins of this weak shock front more particularly. 



3.1. Lack of Supersonic Substructure 

A2556 has a cD galaxy which is consistent with the X-ray peak within 1" as shown in 
Figure [T|d. This indicates A2556 has no violent disturbance caused by major merger in recent 
~ 10^ yrs. The difference in photometric magnitudes between two most luminous galaxies, 
namely the magnitude (luminosity) gap A'magi2, is the simplest measurable statistic to 
quantify the dynamical age of a system (Milosavljevic et al. 2006): larger Amagu indicates 
an older dynamical age. According to the absolute total magnitudes of the brightest (2MASX 
J23130142-2138039) and second-brightest (2MASX J 23132498-2137390) galaxies from Smith 
et al (2004), Amagi2 = 0.851 in R band. Comparing it with the mean magnitude gap of 730 
clusters (Milosavljevic et al. 2006), {Amagi2) — 0.75, we conclude that A2556 is a dynamical 
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old cluster. 

Krick et al. (2007) measured the flux, surface brightness profile, color and substructure 
in the diffuse intracluster light (ICL) of A2556. They concluded that this cluster is dynamical 
relaxed and shows no subcluster. We further examine the line-of-sight velocity distribution of 
46 member galaxies in A2556 identified by Smith et al. (2004) and plot it in Figure HI which 
shows a single velocity plateau from eye's impression. We can fit the velocity distribution 
with a single Gaussian profile, the best fit {x^/dof = 27.87/17) gives the average velocity of 
(v) = 26165.29 ±88.34 km s~^ and the corresponding variance of a = 658.31 ±88.37 km s~^. 
The Kolmogorov-Smirnov statistic tells observed distribution has a probability of 95.0% to 
be a single Gaussian profile, which indicates that A2556 is very unlikely to have subcluster. 

Previously observed shock fronts in ICM of galaxy clusters, for instance bullet cluster 
(Markevitch et al. 2002) and A520 (Markevitch et al. 2005), could be attributed to su- 
personic infalling subclusters. However, there are indications that A2556 possibly has no 
subcluster. Therefore, the weak shock front in A2556 may not be formed by the supersonic 
infalling of the subcluster. 

3.2. Nonlinear Steepening Weak Shock in A2556 

Besides the supersonic flow, the nonlinear steepening from normal acoustic wave could 
also generate shock front (Landau & Lifshitz 1959). Such nonlinear steepening is a possible 
reason for the weak shock front in A2556, since the supersonic infalling of the subcluster 
may not be the origin of the shock front in it. 

Firstly, we would like to discuss the origin of normal acoustic waves in clusters. Accord- 
ing to their numerical simulations Nagai et al. (2003) showed that even for the relatively 
relaxed cluster gas turbulent velocity in the ICM can reach up to about 20% - 30% of the 
local sound speed, and this result is not rule out by current observational upper limits (e.g., 
Sandsers 2011). Such turbulence generates acoustic wave in the ICM (Fujita et al. 2003), 
which can be described by Euler's equation. Since the turbulent velocity v is comparable 
with the sound velocity Cg, the nonlinear term (v ■ grad)v in Euler's equation cannot be 
neglected (Landau & Lifshitz 1959). This leads the acoustic wave to inevitably steepen into 
the shock wave. According to this mechanism, the ICM turbulence of A2556 could be the 
origin of the shock front we detected. 

Secondly, these shock waves would dissipate the energy and heat the surrounding gas 
of the clusters, and we perform a schematic calculation for the nonlinear steepening weak 
shock. Following the method of Stein and Schwartz (1972). When acoustic wave steepens. 
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the crest overtakes the trough in a distance given by 




dL 



X 



(4) 



■s 



2 



where 7 = 5/3 is the adiabatic index for monatomic gas, A is the pulse width, v is the 
pulse velocity amplitude, Cs is the sound speed (for A2556, Cs = 900 km s~^ according to its 
average temperature 3.2 keV), and L is the distance, which the excited sound wave travels 
before it steepens into the shock front. Since energy flux of acoustic wave F oc pv"^ and 
assuming acoustic wave prior to shock formation does not dissipate any energy, we get the 
relation v oc p~^^^. Additionally, with a single (3 distribution p = po[l + {L/rcY]~^^^^ for the 
gas density in the ICM, where Tc is the core radius, we obtain 



where Vq is the initial velocity amplitude. Suppose the acoustic wave originates from the 
cluster center. The distance it travels before steepens to the shock front is L ~ IGO/if/ kpc 
in A2556. According to Eq. [S], A ~ AOOh^i kpc, which means the period of acoustic wave 
generated by the turbulence of ICM in cluster center is ^ 0.43 Gyr. 



The Chandra observation of A2556 reveals a weak shock front located at approximately 
160/2.7/ kpc northeast of the cluster center with Ai = 1.25]1^qq3. Since no subcluster associated 
with the weak shock front is observed in the optical or X-ray band, we cannot ascribe the 
formation of the shock front to the supersonic infalling of the subcluster. We suggest the 
steepening of the acoustic wave could be a possible origin of such shock front, and such 
acoustic wave could induce by the turbulence of the ICM in the cluster core. 

The nonlinear steepening shock front should be an ubiquity phenomenon in the ICM. 
In contrast, the observation of nonlinear steepening weak shock front is rare. This may 
ascribe to that the weak shock wave is difficult to observe. However, to solve this ques- 
tion requires further investigation, including both deeper X-ray observation and extensive 
theoretic studies. 




(5) 



4. 



SUMMARY 
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Table 1. Four Spectra Model Fittings for A2556 in § l2X2l 



Model'' 


T 
(keV) 


Z 

Zq 


T 

(koV) 


M 

(Me yr -1) 


X^/dof 


ITA 

ITM 

FC 

VC 


3.21±0.09 
3.25±0.01 
3.46±0.15 
3.46±0.15 


0.39±0.05 
0.44±0.06 
0.43±0.06 


0.10 


19.99±10.14 
20.13li°3f 


1.13 
1.13 
1.10 
1.10 


''ITA 


represents 


WABSxAPEC model 


; ITM 


represents 



WABSxMEKAL model; For the MEKAL+MCKFLOW model, FC 
represents the lower bound on the gas temperature in the MCKFLOW 
fixed at 0.1 keV and VC represents this component set free. 
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Fig. 1. — (a): Raw Chandra image of Abell 2556 in the 0.7-7.0 keV energy band. Two 
white hnes highhght the range of edge span, (b): SDSS B-band image for the same sky field, 
where the X-ray contours are overlaid. The A2556 image is divided into 4 regions noted as 
E, S, W and N, where SBPs of 4 different directions are extracted. In region E, region 1 and 
region 2 denote the post-shock and pre-shock regions, respectively. White elliptical annuli 
in region E demonstrate regions for extracting spectra in §2.3.1. 
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Fig. 2. — Exposure-corrected SBPs extracted from a series of narrow elliptical annuli across 
the edge of arm-like structure in region E (see Figure lb) and other three regions (S, W and 
N) for comparison. The cross, star cross, circle cross and box cross presents SBPs in regions 
E, S, W and N, respectively. The solid line demonstrates the best fit of the SEP in region 
E according to Eq. [1] A substantial discontinuity of the SEP in region E is located around 
-R=160/if/ kpc; in other three regions, no such discontinuity is observed. 
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Fig. 3. — (a): Deprojected gas temperature profile across the shock in region E. R represents 
the distance from the X-ray peak of A2556. A substantial high temperature is observed 
around i?=160/if/ kpc, which exactly corresponds to the shock front in A2556. (b): The 2- 
dimensional projected temperature map derived from Gu et al. (2009). The white elliptical 
annulus corresponds to region 1 in Figure [TJj. 
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Fig. 4. — Recessional velocity histogram of 46 member galaxies in A2556, which is featured 
with a single Gaussian distribution. The best Gaussian fit gives /dof = 27.87/17 and the 
Kolmogorov-Smirnov statistic tells observed distribution has a probability of 95% to be a 
single Gaussian profile, which indicates that A2556 is very unlikely to have subcluster. 



